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A new single-particle shell model is derived by solving the Schrodinger equation for a semi- 
spheroidal potential well. Only the negative parity states of the Z{z) component of the wave function 
are allowed, so that new magic numbers are obtained for oblate semi-spheroids, semi-sphere and 
prolate semi-spheroids. The semi-spherical magic numbers are identical with those obtained at the 
oblate spheroidal superdeformed shape: 2, 6, 14, 26, 44, 68, 100, 140, ... The superdeformed prolate 
magic numbers of the semi-spheroidal shape are identical with those obtained at the spherical shape 
of the spheroidal harmonic oscillator: 2, 8, 20, 40, 70, 112, 168 ... 

PACS numbers: 03.65.Ge, 21.10.Pc, 31.10.-l-z, 



The spheroidal harmonic oscillator have been used in 
various branches of Physics. Of particular interest was 
the famous single-particle Nilsson model ll| very success- 
ful in Nuclear Physics and its variants for atomic 
clusters. Major spherical-shells N = 2,8,20,40,58,92 
have been found [3] in the mass spectra of sodium clus- 
ters of TV atoms per cluster, and the Clemenger's shell 
model Q was able to explain this sequence of spherical 
magic numbers. 

In the present paper we would like to write explicitly 
the analytical relationships for the energy levels of the 
spheroidal harmonic oscillator and to derive the corre- 
sponding solutions for a semi-spheroidal harmonic oscil- 
lator which may be useful to study atomic cluster de- 
posited on planar surfaces. 

For spheroidal equipotential surfaces, generated by a 
potential with cylindrical symmetry the states of the va- 
lence electrons were found [3| by using an effective single- 
particle Hamiltonian with a potential 
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In order to get analytical solutions we shall neglect an 
additional term proportional to — (l^)n). We plan to 
include in the future such a term which needs a numerical 
solution. K. L. Clemenger introduced the deformation 6 
by expressing the dimensionless two semiaxcs (in units 
of the radius of a sphere with the same volume, Rq = 
rgTV-'^'''^, where is the Wigner-Seitz radius, 2.117 A for 
Na as 
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The spheroid surface equation in dimensionless cylindri- 
cal coordinates p and z is given by 
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where a is the minor (major) semiaxis for prolate (oblate) 
spheroid and c is the major (minor) semiaxis for prolate 
(oblate) spheroid. Volume conservation leads to a^c ~ 1. 

One can separate the variables in the Schrodinger 
equation, H'^ = iSvP, written in cylindrical coordinates. 
As a result the wave function 0, [H may be written as 
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where each component of the wave function is ortonor- 
malized leading to 
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in which r/ = R^p /a\ and the quantum numbers m — 
(rij^ — 2i) with i = 0, 1, ... up to {n±^ — l)/2 for an odd n±^ 
or to {n±^ — 2)/2 for an even n±. L™(a;) is the associated 
Laguerre polynomial and the constant aj^ = yJWJJjuJZ 
has the dimension of a length. 
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where ^ = R^z/az, ctz = \/^/ MlUz , and the main quan- 
tum number n = n±_ + iiz — 0, 1, 2, .... 
The eigenvalues are 
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The parity of the Hermite polynomials _ff„^ (^) is given 
by (—1)"* meaning that the even order Hermite poly- 
nomials are even functions H2nA~C) = ^2nj(C) and 
the odd order Hermite polynomials are odd functions 
H2n^+i{—C) = —H2n-+i{C)- There is a recurrence rela- 
tionship 2zHn = Hn+i+2nHn-i. One has Hq = 1, Hi = 
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FIG. 1: LEFT: Spheroidal harmonic oscillator energy levels in units of Juvo vs. the deformation parameter S. Only 6 major 
shells (A^ — 0, 1, 2, 5) have been considered. Each level is labeled by n, n± quantum numbers and is (2nx +2)-fold degenerate. 
The labels are 0,0; 1,0, 1,1; 2,0, 2,1, 2,2; 3,0, 3,1, 3,2, 3,3; 4,0, 4,1, 4,2, 4,3, 4,4, etc. RIGHT: Semi-spheroidal harmonic 
oscillator energy levels in units of fuvo vs. the deformation coordinate S. Only 9 major shells (A'^ = 0,1,2,..., 8) have been 
considered. Each level is labeled by n, n± quantum numbers (with — n — n± = 1, 3, 5, ... and is {2n± + 2)-fold degenerate. 
The labels are 1,0; 2,1; 3,2, 3,0; 4,3, 4,1; 5,4, 5,2, 5,0; 6,5, 6,3, 6,1, etc. The semi-spherical magic numbers are identical 
with those obtained at the oblate spheroidal superdeformed shape {S — —2/3): 2, 6, 14, 26, 44, 68, 100, 140, ... 



2z, H2 = 4:z^ -2,H3 = 8z^ - I2z, Hi 
H5 = 32z5 - 160z3 -t- 120^, etc. 
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In units of tiioo the eigenvalues, e = E/{Tiu)q)^ are given 
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For a prolate spheroid, (5 > 0, at nj^ = the energy level 
decreases with deformation except for n = 0, but when 
n_L = n it increases. For a given prolate deformation and 
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FIG. 2: LEFT: Harmonic oscillator potential V = l^(C), the 
wave functions Zn^ = Zn^iS,) for Uz — 0,1,2,3,4,5 and the 
corresponding contributions to the total energy levels „^ = 
En^/hiJz = (uz + 1/2) for spherical shapes, 5 = 0. ^ = 
zR^l^hlMuz. RIGHT: The similar functions for a semi- 
spherical harmonic oscillator potential. Only negative parity 
states are retained which are vanishing at ^ = where the 
potential wall is infinitely high . 



a maximum energy e,„, there are firam closed shells and 
other levels for high-order shells up to rimax'- 
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and similar formulae for oblate deformations. (5 < 0. The 
low lying energy levels for the six shells (main quantum 
number n = 0, 1, 2, 3, 4, 5) can be seen in figure [TJ Each 
level, labelled by n± , n, may accomodate 2n± -j- 2 parti- 
cles. One has 2 X]nx=o("-L + 1) = l)(n-|-2) nucleons 
in a completely filled shell charcterized by n, and the 
total number of states of the low-lying n -f 1 shells is 
T,n=oiri+l){n + 2) = {n+l){n + 2){n + 3)/3 leading to 
the magic numbers 2, 8, 20, 40, 70, 112, 168... for a spheri- 
cal shape. Besides the important degeneracy at a spher- 
ical shape {S = 0), one also have degeneracies at some 
superdeformed shapes, e.g. for prolate shapes at the ra- 
tio c/a = (2 + 5)/ (2 - 5) = 2 i.e. 6 = 2/3. More details 
may be found in the Table HI The first five shells can 
reproduce the experimental magic numbers mentioned 
above; in order to describe the other shells Clemenger 
introduced the term proportional to (P — (l^)n)- 

Let us consider a particular shape (half of an oblate or 
prolate spheroid) of a semi-spheroidal cluster deposited 
on a surface with the z axis perpendicular on the sur- 
face and the p axis in the surface plane. Then the semi- 
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FIG. 3: Variation of shell corrections with for Na clusters. TOP: 5 = 0. The semi-spherical magic numbers are identical 
with those obtained at the oblate spheroidal superdeformed shape: 2, 6, 14, 26, 44, 68, 100, 140, ... For a prolate superdeformed 
(S = 2/3) shape the magic numbers are identical with those obtained at the spherical shape: 2, 8, 20, 40, 70, 112, 168, ... 
Other magic numbers are given in Table H] Oblate and prolate shapes are considered on the left-hand side and right-hand side, 
respectively. 



spheroidal surface equation is given by 
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The radius of the semi-sphere obtained for the defor- 
mation 5 = is given by the volume conservation, 
(l/2)(47ri?3/3) = 47ri?i]/3, leading to R, = We 
shall give p, z, a, c in imits of Rg instead of i?o. Accord- 
ing to the volume conservation. a^ci?^/2 = i?g so that 
a?c = 1. Other kind of shapes obtained from a spheroid 
by removing less or more than its half (as in the liquid 
drop calculations 0) wiU be considered in the future; in 
this case it is not possible to obtain analytical solutions. 

The new potential well we have to consider in order to 
solve the quantum mechanical problem is shown in the 
right-hand side of the figure [21 The potential along the 
symmetry axis, Vz{z), has a wall of an infinitely large 
height at z = 0, and concerns only positive values of z 
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In this case the wave functions should vanish in the 
origin, where the potential wall is infinitely high, so 
that only negative parity Hermite polynomials (n^ odd) 
should be taken into consideration. 



From the energy levels given in figure [T] we have to 
select only those corresponding to this condition. In this 
way the former lowest level with ti = 0, nj_ = should 
be excluded. From the two leveles with ?i = 1 we can 
retain the level with nj_ = i.e. = 1. This will be 
the lowest level for the semi-spherical harmonic oscillator 
and will accomodate 2n± + 2 = 2 atoms. From the three 
levels with n = 2 only the one with = nj_ = 1 with 
2n± + 2 — A degeneracy is retained so that the first two 
magic numbers at spherical shape (<5 = 0) are now 2 
followed by 6, etc. Some deformed magic numbers may 
be found in the Table U and as position of minima in 
Fig. El 

Each level, labelled by n^, n, may accomodate 2nj^ -f 2 
particles. When n is an odd number, one should only 
have even n_L in order to select the odd = n — n±. 
The contribution of the shells with odd n to the semi- 
spherical magic numbers will be 
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leading to the sequence 2, 8, 18 for n = 1,3,5. The con- 
tribution of the shells with even n to the semi-spherical 
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TABLE I: TOP: Deformed magic numbers of the spheroidal harmonic oscillator. 
BOTTOM: Deformed magic numbers of the semi-spheroidal harmonic oscillator. 
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2, 8, 18, 20, 34, 38, 58, 64, 92, 
100, 136, 148, ... 
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2, 8, 20, 22, 42, 46, 76, 82, 124, 
134 ... 
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0.4 
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2, 8, 10, 22, 26, 46, 54, 66, 84, 
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2, 4, 10, 16, 28, 40, 60, 80, 110, 
140, ... 


-1 


3 
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4, 12, 18, 24, 36, 48, 60, 80, 100, 
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2, 6, 12, 22, 26, 36, 42, 56, 64, 
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90, 114, 132, ... 


-0.4 
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158, ... 


2/3 
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2, 8, 20, 40, 70, 112, 168, ... 


-1 
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2, 6, 20, 30, 42, 58, 78, 102, 130, 
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2, 8, 10, 14, 22, 26, 46, 54, 66, 
84, 96, 114, 138, 156, ... 
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which gives the sequence 4, 12, 24 for n = 2,4,6. This 
should be interlaced with the preceding one so that the 
magic numbers will be 2, 2 + 4 = 6, 6 + 8 = 14, 14 + 12 = 
26, 26 + 18 = 44, 44+24 = 68, as shown at the right-hand 
side of the Fig. [T] 

The equation ^ from the harmonic oscillator, in units 
of hujo is still valid, but one should only allow the values 
of n and n± for which = n — n± > 1 are odd numbers. 

The ortonormalization condition of the component 
of the wave function became 
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with Uz = 1, 3, 5, n for odd n and Uz = 1, 3, 5, ti — 1 
for even n. Consequently the normalization factor is \/2 
times the preceding one 
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For a nucleus with mass number A the shell ga; 
given by Hivq = 41^1^/'^ MeV. For an atomic cluster 
the single-particle shell gap is given by 
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which is 3.0613iV-i/3 eV in case of Na clusters. Since we 
consider solely monovalent elements, N in this eq. is the 



number of atoms and t denotes the electronic spillout for 
the neutral cluster according to [l3| ■ 



The shell correction energy, 6U Uj , in figure [3] shows 
minima at the oblate and prolate magic numbers given 
in the lower part of the table [D The striking result is 
that the superdeformed prolate magic numbers of the 
semi- spheroidal shape are identical with those obtained 
at the spherical shape of the spheroidal harmonic oscil- 
lator. We expect that this kind of symmetry will not 
be present anylonger for the Hamiltonian including the 
term proportional to (F— (l^)n) and/or the more complex 
equipotential surface we shall study in the future. 
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